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ABSTRACT 


A limiting factor to the development of new suppression technology is the lack of 
appropriate screening methods. Among the new technologies that are intended to replace 
Halon1301'are Solid Propellant Gas Generators, (SPGG’s), or flame suppressing gas 
generators. SPGG’s are a spin-off from airbag technology and have demonstrated their 
ability to suppress certain types of fire, particularly aircraft engine nacelle and dry bay 
fires. To date there is no adequate screening method for SPGG’s. The flow originating 


from burning a solid propellant reduces the Damkohler number by decreasing the 


residence time (high velocity flow of products) and increasing the chemical time (by 
directly altering the reactant concentrations, oxygen displacement effect). If the 
DamkG6hler number decreases below a critical value sudden extinction of the flame 
occurs. A facility that will serve to assess the performance of SPGG’s has to be able to 
evaluate the combined effects of the gas discharge. Several alternatives have been 
proposed for an adequate flame and enclosure that will represent a “worst case” scenario 
for extinction and subsequent re-ignition. Among these alternatives is the recirculation 
zone induced by either a bluff body, a baffle, a backward facing step or a trench inside a 
wind tunnel. These configurations provide a controlled increase in the residence time, 
thus an increase in the Damkéhler number. This report describes a preliminary 
evaluation of these configurations by means of a literature search and some preliminary 
computations using Large Eddy Simulation (LES) code developed at NIST. The 
literature review showed that the proposed configuration is adequate for the present 
application since it creates a recirculation zone with enhanced mixing, entrainment from 
the main stream and product evacuation. The flow parameters can be adjusted to provide 
a stable recirculation zone in a broad range of main stream velocities. Limitations to this 
approach together with a series of design criteria are also proposed. Preliminary 
calculations with the LES code showed qualitative agreement with reported experimental 
data while the literature showed that k-e codes seem inappropriate to model the 
recirculation zone. Different diagnostic techniques used for similar experiments are 
evaluated as candidates for characterization of the experimental facility, emphasis is 
given to velocimetry. Finally, an intermediate scale experimental facility at the LCD- 
University of Poitiers, France, is presented. This fixture has the potential to serve as an 
intermediate scale validation of the proposed screening procedure. 
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1. Background 


When a fire breaks out there is mainly one immediate concern: how to prevent the 
fire from spreading without casualties and with a minimum property loss. For many 
years, at least until 1987 and the signature of the Montreal Protocol, other issues, such as 
the interaction of the suppression process with the local and global environment were a 
matter of little concern. A good example of highly effective suppression agents affected 
by the Montreal Protocol is Halon. Recent environmental problems associated with 
Halons have been well publicized and an amendment to the Protocol required commercial 
Halon production to cease at the beginning of 1994. However, Halons have been the 
agents of choice for numerous fire protection applications mainly because of their great 
ability to inhibit flames at low concentrations but also because of other positive 
characteristics (high liquid density, high stability, low boiling point, toxicity, cost, etc....). 
Consequently large users of Halon have been forced to search for suitable alternatives. 

Following the Montreal Protocol there was a renewed interest in defining new fire 
suppression technology and as a result an early review of the technology available to 
replace Halon was presented in 1993 [1]. DiNenno [2] and Tapscott [3] focused their 
analysis on the search for alternative chemical agents. Grosshandler and al. [4,5] 
discussed research conducted for aircraft applications. New concepts involving chemical 
agents, water-based options, chemically generated gases and aerosols, and physically- 
acting suppression were presented. A summary of the results obtained from extensive 
laboratory-scale experiments was also included. 

It has been experimentally demonstrated that solid propellant gas generators 
(SPGG’s) can be successfully used in a manner similar to a streaming agent to suppress 
fires while remaining environmentally friendly. The combustion products emerging from 
gas generators are considered to have no ozone depletion and global warming potential 
and its use leads to enormous practical advantages such as small size and long storage 
life. Efforts to make the use of this technology possible have significantly increased in 
the last few years with significant involvement of both private and public sectors. 
Notable among these efforts are those of the National Institute of Standards and 
Technology (Building and Fire Research Laboratories), the Navy, Airforce, aircraft 


manufacturers and military contractors. Extensive experimental work has been done by 
the Naval Air Warfare Center Weapons Division in China Lake. In cooperation with 
several aircraft manufacturers, experiments were conducted on real scale facilities (High 
Velocity Airflow System (HIVAS), F/A-18 and V-22 wing dry bay simulators, F/A-18 
engine nacelle simulator) leading to the development of gas generators specifically 
designed for fire suppression (PFE, FSGG). Similar tests have been conducted by the Air 
Force on CFM-56 engines (found on the KC135-R). Details of the status of these 
research programs and other research efforts have been summarized in the proceedings of 
the NIST-1995 Workshop on Solid Propellant Gas Generators [6]. 

Grosshandler et al. [7] pointed that one of the main limitations to the development 
of new suppression technology was the lack of appropriate screening methods. They 
examine two new concepts for testing liquid aerosol and solid propellant gas generator 
(SPGG’s) fire suppression technologies. The first concept is for a bench-scale 
suppression screen suitable to compare the ability of dispersed fluids with differing 
chemical and physical properties to extinguish a laboratory flame; the second is for a 
facility to test SPGG-based agents and release mechanisms. The present review lies in 
the framework of the development of this second concept. 

Definition of a facility that will serve to test SPGG-based agents and release 
mechanisms requires the definition of an adequate flame that will serve as a standard and 
characterization of its interaction with the flow induced by a gas generator. The existing 
knowledge on the characteristics of the flow field in the fire enclosure and the flow 
induced by a gas generator is so far limited to the needs of the current applications 
(cooling air, airbag inflators). In general, the available literature is of global nature and 
consists mainly on average temperature, velocity, pressure and product concentration 
measurements. This information is important as it relates to the nature of the resulting 
atmosphere but provide little insight on the mechanism that lead to flame extinction. 

Since extinction is most likely to occur during the transient stages and re-ignition 
might be possible, before steady conditions are attained. To quantify the potential for re- 
ignition it is necessary to characterize the flow field in the facility and its transient 


behavior after discharge of the SPGG. 


Several alternatives have been proposed for an adequate flame and enclosure that 
will represent a “worst case” scenario for extinction and subsequent re-ignition. It was 
determined by Hamins [8] that a controlled increase in the residence time is obtained 
inside the recirculation zone induced by either a bluff body, a baffle, a backward facing 
step or a trench inside a wind tunnel. Preliminary computations were conducted by 
means of a Large Eddy Simulation (LES) CFD code [9,10] with promising results. The 
base geometry of a wind tunnel with an induced recirculation zone will be explored in 
detail throughout this report and the possibility of using the LES code in the design 
process will be assessed. 

This work is part of a multi-step study including literature survey, experimental, 
theoretical and numerical development of the various aspects of the problem: fluid 
mechanic, combustion and their interactions. This report summarizes a literature review 
and evaluation of a numerical technique conducted with the objective of providing some 
background for the design of a well characterized test fixture for screening the fire 
suppression effectiveness of agents. The facility should be designed with inherent 
flexibility, allowing variation of key geometric parameters that impact the characteristic 
residence time of the recirculation zone (the Damkhdler number). Thus, the literature 
survey focuses primarily on the work related to the description of the structure of the flow 
field behind a backward-facing step or a baffle. Furthermore, a medium scale fixture 
used at the CNRS Poitiers to study the interaction between a cross-flow and a pool fire is 


presented and its ability to serve as an intermediate test facility for SPGG’s is underlined. 


2. Introduction 


Gas generators discharge the cooled down products of a solid propellant 
combustion reaction into ambiance. A gas generator typically consists of a solid 
propellant tablet that will, upon ignition, rapidly react to generate gas phase combustion 
products and particulate, an igniter to initiate the combustion of the propellant, a filter 
system and an exhaust mechanism. Particulate are trapped by a filter that also serves to 
reduce the temperature of the gas-phase combustion products. This element has been 
commonly used as a mechanism to inflate air-bags. The principal gas-phase product of 
the combustion process is nitrogen, therefore, a logical extension to this technology is fire 
suppression [6]. 

If a flame is subject to a fast flow of nitrogen extinction might occur. The 
advantages of generating this fast flow by means of a solid propellant are many: it is 
compact, it has a very long storage and service life, it can be used in areas of difficult 
access, the response time can be extremely fast and this type of system is considered to 
have no ozone depletion or global warming potential. One of the main issues that has 
slowed the development of gas generators as a fire suppression technique is the lack of 
adequate test protocol to assess their performance. Estimation of the reliability and 
efficiency of this technique depends on a better understanding of the high speed flow and 
the fundamental chemical and thermal processes involved in extinction and re-ignition of 
a fire. 

The process of extinction can be described by means of the Damk6éhler number 
(Da=(Residence Time)/(Chemical Time)), by either reducing the residence time or 
increasing the chemical time a critical Damkohler number for extinction can be attained. 
The flow originating from burning a solid propellant can reduce the residence time (high 
velocity flow of products) and increase the chemical time (by directly altering the reactant 
concentrations, oxygen displacement effect) resulting in sudden extinction of the flame. 


Since the changes induced by the burning solid propellants are very strong and sudden, 


the extinction process is expected to be extremely fast. This suppression technique has 


potential applications as a substitute to Halon 1301 in engine nacelles, dry bays and army 
vehicles as well as localized and difficult to access fires. 

Extinction mechanisms have been a subject of numerous studies and excellent 
reviews can be found in the literature. A. good summary of the existing knowledge is 


provided by Williams [11]. The process of extinction can be described by means of the 


Damkoéhler number (Da=(Residence Time)/(Chemical Time)), by either reducing the 


residence time or increasing the chemical time a critical Damkoéhler number for extinction 


can be attained. The flow originating from burning a solid propellant can reduce the 
residence time (high velocity flow of products) and increase the chemical time (by 
directly altering the reactant concentrations, oxygen displacement effect or by decreasing 
the temperature of the pyrolysis and reaction zones) resulting in sudden extinction of the 
flame. Furthermore, the interaction between the flame and the flow, coming from the gas 
generator, can affect the turbulent structure of the system. Local stretch rates can be 
altered by the nitrogen rich flow resulting in local extinction zones. Depending on the 
geometry of the system, the type of flame, the characteristics of the gas generator and the 
distance between the gas generator and the flame the effects of each independent 
mechanism can be very different. 

An issue of great relevance to this specific problem is the relative importance in 
the extinction process of radical depletion and temperature decrease. The flow induced 
by the gas generator will affect the chemistry of the reaction as well as decrease the 
temperature of the flame and fuel. Although the imposition of a non-oxidizing flow and 
the enhanced heat losses induced by the gas generator tend to, both, favor the extinction 
process, their relative importance will have a strong effect on the range of scenarios 
where a gas generator can be used as a fire suppression mechanism. Geometrical 
considerations will be different for both extinction mechanisms. Significant work on the 
interaction between flames and high speed flow has been reported in the past and have 
been reviewed by Williams [11] and Blazowski [12] but all the information pertains to 


scenarios that significantly differ from that of fire flame extinction. 


When considering the sudden onset of an enhanced velocity field past a given 
object it can be observed that the transient conditions under which the object will 
continue to burn (or extinguish) differ little from those under which it will continue to 
burn (or extinguish) in a steady field of equal intensity and distribution. If the flame can 
not persist after the sudden onset of the velocity it is likely that it will not re-ignite once 
steady conditions are attained. Cooling down of the fuel will then follow [12]. In the 
case of solid propellant induced flows the extent and characteristics of the transient and 
steady state conditions depend on the geometrical configuration. The relative position of 
the fire with respect to the generator will significantly affect the flow structure and 
development. If the fire is close to the generator and the characteristic volume of the 
geometrical boundaries are very large compared to that of the fire and the gas flow 
induced by the generator, then both the transient and steady state periods will be very 
short and the flow field will regain its initial characteristics very fast. If the characteristic 
volume of the flow created by the gas generator is comparable to the volume of the room 
then a short transient period will be followed by a long steady period (flooding), in this 
case the flow field will not be expected to regain its initial characteristics. Obstructions 
in the path between the generator and the fire will significantly affect the velocity and 
distribution of the flow reaching the fire and will also affect the extent and nature of the 


transient process. 
If flooding occurs it is expected that the critical Damkohler number will be 


attained permanently and cooling down of the fuel, which is a very slow process (when 
compared to the chemical and residence times of the flame), will follow. If both transient 


and steady periods are very short, first scenario depicted above, the flow conditions will 
provide a DamkGéhler number that will allow re-ignition if enough gaseous fuel can still 


be produced when the low velocities and high oxygen concentrations are re-instated 
[13,14]. Fuel pyrolysis will depend on the surface temperature and for many materials, 
when burning, the surface temperature is much higher than the vaporization temperature 
(i.e. charring materials), therefore, fuel pyrolysis will occur even after extinction of the 


flame. Comparison between the characteristic time of the flow induced by the gas 


generator and the characteristic cooling time becomes relevant and re-ignition becomes a 
significant issue when evaluating the performance of gas generators as a suppression 
mechanisms. 

A simple way to increase the residence time is by creating a recirculation zone, 
this will significantly reduced the strain on the flame and thus provide a “worst case 
scenario” to test gas generators. The different issues that need to be addressed when 
designing such a facility will be covered throughout a review of relevant studies. These 
studies include work on the flow structure behind a backward facing step, a baffle and a 
trench. The review is complemented by some background information on gas generators 
and preliminary CFD results using the LES [9,10] code and validated against 


experimental data obtained from the literature survey. 


3. Solid Propellant Gas Generators (SPGG’s) 


Solid propellant gas generators, (SPGG’s), or flame suppressing gas generators, a 
spin-off from airbag technologies, have demonstrated their ability to suppress certain 
types of fire, particularly aircraft engine nacelle and dry bay fires [15]. As stated by 
Yang and Grosshandler [6] the main areas where research will contribute to the efficient 


implementation of this technology are: 


e Identification of certification procedures for gas generators in fire suppression 
applications, 

e Determination of critical parameters for evaluating fire suppression applications. 

e Development of a standard methodology to facilitate testing of gas generators. 

e Identification of possible applications other than protection of engine nacelles and dry 
bays, 


e Identification of a new generation of propellants. 


According to Yang and Grosshandler [6], a SPGG is essentially an airbag inflator 
without a bag. That is, the gas generated is discharged directly into ambience rather into 


a bag. A typical SPGG consists of: 


e a solid propellant tablets which will, upon ignition, rapidly react to generate gas- 
phase combustion products and particulate 

e an igniter to initiate the combustion of the propellant 

e a filter system to prevent or minimize the release of particulate from the combustion 
reactions into the ambience 

e a heat transfer system to cool the high temperature combustion gas before being 
discharged into the ambience 


e an exhaust mechanism to disperse the gas efficiently. 


There are basically two types of airbag inflator systems: the conventional and the 
pre-pressurized or gas-assisted. In a conventional system, the gas that is used to inflate 
the bag depends entirely on the combustion gas generated by the solid propellant. 
However, in a pre-pressurized or gas-assisted system, the high temperature gas as a result 
of the combustion of the propellant is first mixed with a pre-pressurized inert gas at 
ambient temperature before being discharge into a bag. Similarly, one can also 
conveniently classify solid propellant gas generators into two categories, depending upon 
their functions: conventional and hybrid. When a gas generator is used alone for fire 
suppression, it is termed “conventional”. When it is used together with other liquid or 
powdered fire suppressing agents, it is termed “hybrid”. In a hybrid system, the gas 
generator normally is used as a means to provide sufficient pressurization so that the 
expulsion of liquid or powdered agent from a storage vessel can be facilitated. A 
different alternative or hybrid inflator is composed of pressurized argon that is heated by 
a small charge to produce a sufficient volume of gas for bag inflation [16] 

The attractiveness of using SPGG’s in fire suppression applications lies in the fact 
that the system, when used alone, is considered to have no ozone depletion and global 


warming potential, and is physically very compact. Being a derivative from the airbag 


inflator technology, there are numerous research publications available. Another 
advantage is that since gases are generated via solid propellant reactions, the system can, 
in principle, be tailored to function over a period of few milliseconds to few seconds by 
manipulating the parameters that control the combustion mechanisms. In addition, the 
gas generators have very extended storage and service life. However the toxicity of some 
of the by-products cannot be ignored. 

Sodium azide is the current principal chemical used in solid propellants for gas 
generators. Because of its potential health hazards airbag manufacturers have focused 
recent research on the “non-azide based” propellants. The development of any new 
propellant should include the improvement of: propellant thermochemistry and 
stoichiometry, propellant ignitability and burning rates under various conditions, toxicity 
of combustion products, stability of propellant during storage and transport, and 
propellant thermal properties. In addition issues that need to be addressed are the grain 
size, shape of the propellant and how the propellant is packed in the gas generator. 

Most relevant information related to SPGG’s is present in the airbag inflators and 
solid propellant literature. Global characteristics of the performance of automotive 
airbags address issues specific to this application [16-19] but might be relevant to fire 
suppression. One of these issues is the rate at which the combustion products are 
released. Generally this rate is controlled by cutting the sodium azide propellant into 
pellets. The size of the pellets determines the surface area, which subsequently 
determines the rate of gas generation. Typical values for gas generation rates result in 
airbags being inflated in 20 to 40 milliseconds [19]. Also related to the gas generation is 
the nature of the combustion process. Berger and Butler [20] examine the performance of 
gas-generating propellants by comparing the theoretical combustion behavior of three 
condensed-phase propellants commonly used in the airbag industry. They considered an 
azide and a non-azide propellant, and a double-base propellant. They investigated 
various thermophysical properties including the flame temperature and chemical 
composition of the product gases, the number of gaseous moles produced per mass of 
propellant consumed, the condensed-phase or slag production of each propellant, and the 


toxicity of the combustion products. These investigations were carried out by means of 


the widely used test tank. This work concludes that the thermophysical properties of the 
propellant have a significant effect in the performance of the SPGG. 

Average experimental measurements of the output of a gas generator have always 
been of questionable validity, mainly due to the unsteady nature of the process [20]. 
Many attempts have been made to provide local and transient measurements, among the 
different methods proposed, Wang [21] developed a new semi-analytical procedure, 
called the dual-pressure method, for computing the output from a pyrotechnic inflator. 
This method appeared better than the more traditional average temperature method, since 
it predicts the transient gas temperature in inflators. Although the transient evolution of 
the pressure can be tracked, the results remain of qualitative nature. 

In an attempt to cover the lack of quantitative experimental work modeling has 
been extensively used to evaluate the performance of SPGG’s. Motevalli and Bedewi 
[22] presented a general overview of the modeling effort in the area of airbag 
deployment. This review describes the main possibilities of existing codes and the 
required developments necessary to improve them. This paper points that numerical 
codes use uniform gas thermodynamic properties and empirical data to calculate gas 
generation and calculation of an airbag internal pressure. To obtain an adequate 
characterization of the effluent it is necessary to have the ability to conduct a time- 
dependent calculation of the gas flow into and out of the airbag that includes pressure 
interactions. In addition, gas jet effects need to be incorporated which forces the solution 
of the full Navier-Stokes equations. Furthermore, it was recommended that a sensibility 
analysis of the current models should be conducted. 

Work has been conducted in some of the areas considered of importance by 
Motevalli and Bedewi [22] but the results are still far from being complete. Lupker and 
Bruijs [23] and Groenenboom et al. [24] present modeling results that concentrate on the 
jet issuing from the gas generator and its impact during the airbag inflation. A physical 
model for a conventional gas generator and a gas assisted pyrotechnic inflator was 
developed and validated by Krier and Butler [25] and Butler et al. [26]. 

Butler et al. [26] developed a mathematical model to simulate the transient, 


thermochemical events associated with ignition and combustion of pyrotechnic 


automotive airbag inflator. The conservation equations for mass and energy were derived 
for the interior combustion chamber, filter/cooling screens, exterior plenum, and 
discharge tank. After a brief description of the model development and of the physical 
assumptions made in the analysis, two series of test calculations were presented. The first 
series of calculations is considered as the baseline case of a conventional pyrotechnic 
inflator system that is characteristic of a standard discharge tank validation experiment. 
Transient pressure and temperature profiles generated by the airbag inflator model are 
presented along with properties at the exit nozzles. The parametric study demonstrates 
the usefulness of airbag inflator simulation in assessing the sensitivity of airbag pressure 
curves to various design parameters such as propellant and hardware properties and 
hardware dimensions. The second series of calculations illustrates the influence of pre- 
pressurized inert gas on the performance of a pre-pressurized pyrotechnic inflator system. 
Performance of the inflators is given in terms of pressure-time and temperature-time 
profiles in the inflator and discharge tank as well as pressure-time integrals at specified 
time after ignition. The pre-pressurized pyrotechnic inflator shows certain advantages 
over conventional pyrotechnic units, including significantly lower requirements for solid 
propellant mass, lower operating temperature, more uniform performance at hot and cold 
ambient conditions, and higher thermal efficiency. The chemical composition of the inert 
gas, pre-pressurized system, is also shown to influence the working process of inflator. 
Validation of these results was made with restricted experimental data. Further 
development and experimental validation of this model could result in a very useful tool 
for the design and the development of new inflators, new pyrotechnic compositions, 
ignitors and filters and also of the propellant grain structure. 

In a recent paper Schmitt and al. [27] extend the previous numerical work to study 
the transient operation of a non-azide propellant pre-pressurized airbag inflator. 
Predicted performance of the inflator is presented in terms of pressure, temperature and 
mass flow rate profiles in the inflator and discharge tank which is used to simulate an 
airbag. This work also predicts first-order estimates of gas-phase species exit 
concentrations and characteristic residence times in the inflator. Carbon monoxide, 


produced as a product of combustion from the high flame temperature propellant, is 


partially converted to CO, as it flows from the internal combustion chamber to the 
pressurized plenum before being discharged into the airbag. Specifically, the 
production/destruction of CO is tracked using three different gas-phase reaction models: 
chemically frozen, local equilibrium and finite-rate elementary kinetics. The results 
unambiguously express the need for an airbag combustion program that includes finite- 
rate, gas-phase kinetics. Finally, the results from the finite-rate CO chemistry model 
seem in qualitative agreement with the experimental results reported in the literature. 
Computer codes that are currently used to simulate airbag inflator performance 
address almost exclusively the internal performance of airbag inflators. Chemical 
equilibrium is assumed to determine the products of combustion and flame temperatures. 
Since the gas generation processes are extremely rapid and over in such a short duration, 
chemical equilibrium may not be reached, and simplified or detailed chemical kinetics 
should be considered in the future. In addition, the interaction of the exhaust gas from the 
gas generator with the ambience has to be taken into account. Addressing the above 
mentioned issues could allow to use these computer codes to evaluate the performance of 


SPGG’s. 


4. Literature Review 


As mentioned before, a common way to increase in a controlled way the residence 
time is by creating a recirculation zone. If the recirculation zone is established at the 
flame location a “worst case” fire scenario is effectively created. Recirculation zones can 
be created in many ways and different parameters, including the flame, can affect its size 
and stability. A literature survey was conducted to gather relevant information that will 
lead to the definition of the most effective recirculation zone for the purpose of screening 


SPGG’s. 
a. Structure of the Flow-Field Behind a Step 


A review of previous experimental, numerical and theoretical studies will be 


presented. The flow behind a step corresponds to a typical reattachment process. It is 


12 


worth noting that the reattachment of a shear layer is an important process in a large 
number of practical engineering configurations including, diffusers, buildings, airfoils, 
flame stabilization in engines and combustors. Among two-dimensional flows, the 
backward-facing step is the simplest reattaching flow. The separation line is straight and 
fixed at the edge of the step, and there is only one separation zone instead of two, as seen 
in a flow over an obstacle. In addition, the streamlines are nearly parallel to the wall at 
the separation point, so the step only has a significant effect downstream of the separation 
point. 

The literature presents a large number of studies on large Reynolds number 
turbulent flows but only a few papers dealing with the structure of a laminar separation- 
reattachment region of the flow behind a step. A pioneer paper is due to Acrivos and al. 
[28] who presented results for the steady separated flow past a variety of bluff objects. 
Among the results there is a pressure distribution along the base of a backward-facing 
step. Goldstein et al. [29] conducted an experimental investigation of the laminar flow of 
air over a downstream-facing step. They reported results over a range of 0.36 - 1.02 cm 
in step height and 0.61 - 2.44 m/s in free stream velocity at the step. The main 
observation was that laminar reattachment length is not a constant number of step heights 
as for turbulent flow, but varies with Reynolds number and boundary layer thickness at 
the step. The shape of the velocity profile at reattachment is found to be similar to the 
shape of the velocity profile at separation and downstream of reattachment. 

Bradshaw and Wong [30] provided the first review of the experimental data for 
reattaching flows. The review includes experiments on the low-speed flow downstream of 
steps and fences and especially downstream of a backward facing step. They demonstrate 
the complicated nature of the flow in the reattachment region and its effect on the slow 
non-monotonic return of the shear layer to the ordinary boundary-layer state. A key 
feature of the flow observed by Bradshaw and Wong was the splitting of the shear layer 
at reattachment, where part of the flow is deflected upstream into the re-circulating flow 
region to supply the entrainment. Experiments conducted with low velocity turbulent 
flow [31] showed that at reattachment the re-circulating flow differs from the plane- 


mixing layer in one important aspect: the flow on the entrainment side is highly turbulent, 


as opposed to the low turbulence-level downstream of the reattachment point. The re- 
circulating flow cannot be characterized as a dead air zone. The maximum measured 
back-flow velocity is usually over 20% of the free stream velocity and very large negative 
friction coefficients are observed [32]. The flow in this region is very unsteady with very 
large turbulent structures (at least as large as the step height) passing through the 
reattachment region [33,34]. These observations are of great importance in the present 
application since it ensures the entrainment of air to the re-circulation zone, adequate 
mixing and the evacuation of combustion products. 

The unsteady nature of the flow has been described by means of flow 
visualization which showed that the length of the separation region fluctuated, therefore 
the impingement location of the shear layer moved up and downstream. Quantitative 
measurements confirmed this conclusion and showed that the short-time averaged 


reattachment point deviated from the long-time averaged reattachment location by as 


much as + 2 step heights. A non-dimensional frequency for this motion (f ) was is given 


by the following expression 


where f is the fluctuation frequency X, is the time averaged distance between the 
reattachment location and the step and U, the inviscid-flow velocity upstream of the step. 


The characteristics of the flow seemed to be governed by this non-dimensional parameter. 
Mulin and al. [35] reported that for f<0.5 the recirculation appeared to wash out 


periodically. For f>0.5 the recirculation zone appeared stable but smaller than for the 
steady free-stream case. Further experimental results that tracked the maximum energy 
content of wall pressure fluctuations [36] validated this observation. 

Flow reversal in the reattachment region appears to occur in a random manner 
with no apparent correlation between the near-wall flow upstream and downstream of 
reattachment. The largest structure in the flow originates from the roll-up and multiple 


pairing of span wise vortices [37]. This roll-up is similar to the vortex roll-up and a 
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pairing process is seen in the free-shear layer. The convective speed of these structures is 
about 0.6U,. The span wise coherence or organization of these vortical structures starts to 
break down about 3 step heights downstream of detachment with the turbulence structure 
becoming fully three-dimensional upstream of reattachment. The turbulence intensity of 
the detached flow is 5-10% higher than for the plane-mixing layers. This is believed to 


be the result of a very low frequency (f<0.1) vertical or flapping motion of the 

reattaching shear layer. 

The vertical or flapping motion of the reattaching shear layer encourages the 
believe that this kind of detached flow can be manipulated by imposed unsteadiness. 
Reisenthel and al. [38] used an oscillating flap of length “h” at the backstep. At low 
frequencies (k= f.h/U, = 0.01-0.05) and small flap extensions (45°), the flap does not 
generate much vorticity but periodically blocks the accumulation of vortical fluid at the 
backstep. This causes the shear-layer vortices to be drawn closer to the plate, thus 
shortening the reattachment length. For k= 0.06-0.09 and full flap extension (90°), the 
flap tip produces significant vorticity that merges with the shear layer to reduce X, by 
50%. Beyond k=0.1 with a full flap extension, the flap vortices are so close together that 
they do not allow the separated flow to reattach until they have diffused downstream, 
which produces a larger X,. 

The impingement point of the shear layer moves up and downstream as a function 
of time [32,33,34] but, average values have been commonly used to characterize the 
reattachment zone. The most important dependent parameter characterizing the flow 
field seems to be the average reattachment length which varies from 0 to more than 20 
step heights (Figure 1). Comparison of the reattachment length from various experiments 
provides insight into the effects of varying the following six principal independent 
parameters: 

1. The initial state of the flow. As an example Figure 1 shows the differences between a 
parabolic and a fully developed velocity profile immediately upstream of the step 
(Ref.[41]-p) and a developing velocity profile [Ref. [41] ). As the flow becomes fully 
turbulent, Re,>6000, (Ref. [39]) the length of the recirculation zone becomes 


independent of the Reynolds number and other recirculation patterns appear at the top 


and bottom of the channel [39]. According to Bradshaw [42], recirculation zones 
originating from laminar boundary layers grow more rapidly than those originating 


from turbulent boundary layers. 
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Figure 1. Normalized size of the recirculation zone as a function of the Reynolds number 


. The pressure gradient can have a significant effect on the reattachment length [43]. 

. The aspect ratio of the test section has a negligible effect on the re-attachment length 
for aspect ratios greater than ten [44]. Below ten the reattachment length increases for 
a laminar boundary-layer and decreases for a turbulent flow, the difference is 
apparently due to the effects of the corner flow. Figure 1 shows experimental results 
for an aspect ratio of 1.2 (Ref.[40]) and aspect ratio of 10 (Ref. [39]) showing the 
above mentioned differences. 

. Blowing or suction tunnels behave in a very similar way. 

. Blowing immediately downstream of the step (for flame attachment purposes) was 
studied by Richardson et al [45] in a section with aspect ratios of six and four, flow 
velocities between 50 and 80 m/s and injection velocities ranging from 0 to 0.5 m/s. 


It was observed that the recirculation zone changes in size and injection had a de- 
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stabilizing effect. This results correspond to an entirely different velocity regime but 
the results will be considered when defining the present study. 

6. The Reynolds number. At very low Reynolds numbers, the flow is viscous, and the 
recirculation zone is formed of one stationary eddy at the corner of the step. As the 
Reynolds number increases, an eddy may detach and decay as it moves downstream. 
At moderate values, the flow reaches a state of transition, with eddy shedding at the 
step. At high Reynolds numbers, the flow exhibits turbulent behavior with a 
continuous process of eddy formation and pairing. The length of the recirculation 
increases with Reynolds number, reaching a maximum at transition and then decay to 
a shorter length at the turbulent range. Turbulence statistics show a rise at transition 
(Figure 1). Flow visualization showed that the length of the separation region 
fluctuated so that the impingement location of the shear layer moved up and 


downstream [30]. 


Numerical calculations incorporating the k-e and algebraic-stress model terribly over 
predict experimental observations and especially the turbulence dissipation rate [46]. The 
numerical study of Amano and Goel [47], using the same principles showed similar 
problems. In both cases the recirculating region (size and characteristics) is poorly 
described. Miau and al. [48] indicated that a standard k-e model can yield good 
predictions for the back-step problem only when modified with an independently 
calibrated anisotropic eddy viscosity. Preliminary computations conducted using the 
methodology described by McGrattan et al [9], on a one aspect ratio section, show 
qualitative agreement with the experimental data (Figure 1 [Model]), further exploration 
of the use of this model needs to be done and is part of the present proposal. The results 
presented in Figure 1 do not correspond to the average length of the recirculation zone 
but to the arithmetic average between the maximum and minimum transient values, 
therefore the discrepancy between the model and the experimental results still requires 
further analysis. Quantitative measurements have shown that the short-time average 
reattachment differed from the long-time average reattachment location by as much as + 2 


step heights. 


Many of the papers above include comparative studies where the backward facing 
step was substituted by either a baffle or a trench. Other studies include sudden 
expansions and contractions. The qualitative observations presented above seem to be 
general to all configurations and quantitative details are Acrivos et al [28] Bradshaw and 
Wong [30] and Eaton and Johnston [31]. Emphasis has been given to the backward facing 


step because among two-dimensional flows this is the simplest form of reattaching flow. 


b. Heat Transfer Characteristics Behind a Step 

Heat transfer in the separation and reattachment zones has been the subject of a 
very reduced number of studies [49-58]. Earlier work by Seban [57] experimentally 
obtained local heat-transfer coefficients as a function of the flow velocity. Wall 
temperature distributions in the separated and reattached regions of the flow were also 
presented. It was found that neither region presented similarity patterns common to flows 
that are dominated by the friction at the wall. The effect of suction or injection through a 
slot at the base of the step was also studied indicating small effects on both the pressure 
distribution and the local heat-transfer coefficient. More recent work by Richardson et al. 
[45] has shown contradictory results. 

Filetti and Kays [56] presented experimental data for local heat transfer rates near 
the entrance of a flat duct in which there is an abrupt symmetrical enlargement in flow 
cross section. Although the Reynolds number used for this experiments 70,000<Re,< 
205,000) do not correspond to the present application, the results are of interest in that 
they show that maximum heat transfer occurs at the point of reattachment, followed by a 
decay toward the values for fully developed duct flow. Empirical equations are obtained 
for the Nusselt number at the reattachment point as a function of the Reynolds number. 

The only work in low speed wind tunnel is that of Aung [58] who presented 
experiments conducted with a contraction ratio of 12.5/1. Temperature distributions in 
the gas were measured using a Mach-Zehnder interferometer. Heat transfer upstream of 
the step is shown to be strongly enhanced by streamline curvature. Downstream of the 
step the heat transfer increases monotonically in the streamwise direction but is always 


less than the flat-plate value. Shishov and al. [49] developed an approximate calculation 


method for predicting velocity and temperature profiles, reattachment length and local 
friction and heat transfer coefficients. Their derived analytical expression for heat 
transfer in the reattachment region correlated well experimental heat transfer data for 
different types of flow with a sudden expansion. 

The work of Yogesh and Raghunandan [50] is devoted to the study of the flow 
structure and heat transfer characteristics behind a diaphragm in the presence of a 
diffusion flame. The diffusion flame is stabilized by injecting gaseous fuel into a main 
stream of air or N,-O, mixtures. When a diaphragm obstructs the flow over a surface, the 
flow gets separated at a point upstream of the diaphragm and gets reattached further 
downstream. Consequently, there are two recirculation zones one in the front and the 
other behind the diaphragm. The emphasis of this paper is, mainly, on the nature and 
effects of the downstream recirculation. This attempt to investigate the complex flow 
resulting from combustion of a fuel surface downstream of a diaphragm has led to some 
important results concerning large scale features such as heat transfer rates to the surface. 
Finer features related to the flow pattern and temperature distribution were obtained by 
particle tracking and fine thermocouples. Blowing or fuel injection through the porous 
plate tends to lift and stretch the recirculation zone. The flame is located at the shear 
layer and the eddies appear to be further enlarged. Temperature profiles indicate the 
expected well-stirred type of regions behind the diaphragm with the surface gradient 
temperature gradually increasing along the flow direction. The heat flux to the surface 
shows a constant flux zone near the diaphragm followed by a monotonic rise. The two 
zones can be attributed to the two different eddies visualized by particle track 
photographs. Experiments with a heated main stream with no combustion and the flow 
visualization in the cold flow lead to the conclusion that the flow pattern and heat transfer 
trends are similar with or without combustion reactions. The presence of a diaphragm 
helps in postponing the blow-off of the flame in relation to the boundary layer diffusion 
flame. However, the stable flame zone marked by a perfect two-dimensional flame is 
only marginally widened. In the boundary layer case, due to the proximity of the flame, 
the heat transfer rates are high around the leading edge but sharply decrease along the 


flow direction. 


c. Ignition of a Flammable Mixture Over a Hot Plate 

It is not the objective of this section to summarize the extensive literature on the 
subject of ignition of a flammable mixture but just to present the main parameters that 
need to be considered when studying this process. Most of the work on this subject has 
been conducted by using boundary layer flow and has been summarized in detail by Law 
and Law [59] and Trevino and Sen [60]. Theoretical work has mostly used asymptotic 
analysis (high and low activation energy) to determine the parameters that affect the 


ignition process. The following parameters have been found to affect determine ignition: 


e Critical DamkGéhler number; characteristic residence (flow characteristics) time and 
chemical times (fuel and oxidizer, fuel and oxydizer concentrations, gas pahase 


temperature, etc.) determine a critical Damkohler number that needs to be attained for 


ignition to occur. The attainment of the critical Damkohler number will determine 


the location of the ignition point. The fuel concentration in the flow is highly 
dependent on the pyrolysis rate of the fuel and its nature (charring, non-charring). 

e The thermal properties of the flammable mixture. 

e The thermal properties of the plate; Trevino and Fernandez-Pello [61] show that the 
thermal properties of the material influenced the hot plate influenced ignition. An 
isothermal wall can be considered as the best case scenario for ignition, the thermal 
properties of the wall will indicate how restrictive heat transfer from the wall to the 
flammable mixture is and subsequently how it deviates from the ideal case. 

e Catalytic reactions at the surface; many papers have dealt with the effect of catalytic 
reactions at the hot plate surface. The characteristic chemical time of the catalytic 
reactions is the main parameter that quantifies the importance of the chemical 
interaction between the solid and gas phases. Slow catalytic reactions have negligible 
effect on ignition, fast catalytic reactions have a dominant effect (favorable or 


adverse) on ignition [62]. 
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5. Diagnostic Techniques 


The specific constraints of the present application result in a reduced number of 
diagnostic techniques that have been proven to provide information that could be useful 
in the characterization of the controlling mechanisms of extinction by means of SPGG’s. 
This summary only includes those techniques that have been used in studies related to 
_reattaching flows and concentrates on techniques that provide information on the 
structure and global properties of the flow. Determination of the diagnostic techniques 
that will determine species concentrations would be of significance but goes beyond the 


scope of this review. 


a. Velocity Measurements 


Different diagnostic techniques have been used to characterize the flow field behind a 
backwards facing step. The general criteria for the determination of the adequate 
diagnostic technique has been the Reynolds number. For the low Reynolds number 
regime (Re,<2,000) flow visualization by means of illumination of particles introduced in 
the stream has proven to provide best results. Finaish and al. [63] successfully 
investigated the development of two-dimensional vortex patterns in an accelerating flow 
by means of streak-line visualization. Liu et al. [64] used Particle Image Velocimetry 
(PIV) to study the transient behavior of a turbulent (Re=2,872) flow inside a channel of 
aspect ratio of 12. Bandyopadhyay [65] studied the nature of instabilities and large 
structures generated in reattaching boundary layers and especially behind a backward- 
facing step by the simultaneous use flow visualization and hot-wire anemometry. Plane- 
flow visualization was conducted in the backward-facing step showing that the three- 
dimensional mixing layer developing after detachment persists even after reattachment. 
Velocity measurements obtained by means of both techniques proved to be identical with 
flow visualization having the advantage of providing a flow field instead of a single 
measurement. Honji [66] investigated the incompressible starting flow past a 


downstream-facing right-angled step at Re,<500 by means of flow-visualization 
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techniques. The starting flow down a step height of 1.85 cm was characterized by the 
Reynolds number based on the step height (Re,=U.S/v) and the dimensionless time, 
(t =U.t/S). The formation of secondary vortices, in the recirculation zone, for Re,>120 
and a third vortex for Re,>140 was characterize by means of this technique. In a similar 
way, Sinha and al.[40] reported the results of an experimental investigation carried out at 
a wind speed of 1.8m/s over 3 backward facing steps 0.625, 1.25, and 2.5cm in height. In 
all cases the separating boundary layer was laminar and 1.4cm thick. Data was obtained 
by means of smoke visualization and hot wire measurements. The data served to further 
characterize the reattachment length in the Reynolds number range 100 — 12,500 and to 
describe the general patterns of the flow between the point of separation and that of 
reattachment. 

An important consideration for the success of PIV measurements is determined by 
type of material used for seeding the flow. The constraints placed on the seeding material 
are: 1) the particle’s inertia must be small enough that it faithfully follows the motion of 
the fluid, 2) the particles settling velocity must be small in comparison to the typical 
velocities of the flow, and 3) the particle must be large enough so that it will provide 
sufficient scattered light intensity to register on the CCD [67,68]. The best candidates for 
the current application would be fine hollow glass beads (~ 2 um diameter) or Titanium 
Oxide (TiO,) power (~ 4 pm) introduced through an aerated bed. The suitability of these 
particles for the imaging conditions above can be checked by examining hydrodynamic 
response of the particles to oscillatory forcing, checking the settling velocity, and 
estimating the expected scattered light intensity from the particles. 

The response of a particle to unsteady motion has been calculated by Hjemflet and 
Mockros [69], and for large particle to fluid density ratios (p, /p, <<1) they found the 
expression 
u 1 


Ur /1+(2nSty 
where St = D’p,f /18p,v is the Stokes number and is equivalent to the ratio of the 


particle viscous response time to the frequency of fluctuation of the flow, f The 
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calculated frequency response (where u, /u, = 0.99) for the glass beads and the TiO, is 
then found to be approximately 2300 Hz and 115 Hz, respectively. These values are in 
excess of any unsteady fluid motion expected for the recirculating flows of interest. 
The settling velocity can be estimated using a simple balance between the Stokes 
drag force and buoyancy, giving 
HDIP;S 
pean Sry 


The calculated settling velocities correspond to 0.01 cm/s for the glass beads and 0.2 cm/s 
for the TiO, both of which are small compared to the expected mean velocities of 0.1 — 
1.0 m/s. 

The last quantity of concern is whether the particles scatter enough light to be 
detected by the camera. Adrian and Yao [70] have made calculations of mean exposure 
level generated by spherical particles under different imaging conditions, using the 
equation 
ee Awe) ror [ord 


where A is the wavelength of the light, M is the magnification, d, is the diffraction limited 


point-response diameter of the lens, W =I,6t is the energy of the incident laser pulse, Ay is 
the width of the light sheet, Az is the light sheet thickness, o is the Mie scattering 
coefficient, and Q is the solid angle subtended by the imaging lens. The most stringent 
conditions will occur for the smallest particles at large f#, and using Adrian and Yao’s 
calculations (see their figure 6) give a non-dimensional exposure intensity of 
eAyAz/W=7x10° for D = 2 um and f#=22. For the 2 W laser formed into a 10 
cmx 0.1 cm light sheet used in the proposed experiments, this gives a dimensional 
irradiance of 3.5 nW/cm’, which is 35 times greater than the minimum irradiance 
required for the camera (0.1 nW/cm). Further gains in light intensity can be made by 


lowering the f# several stops. 
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b. Time Dependent Density Measurements 

The use of Mach-Zender Interferometry, Schliren Deflectometry and Rainbow 
Schliren Deflectometry is documented in the literature. The results obtained are generaly 
of qualitative nature. Rainbow Schliren Deflectometry offers significant advantages over 
traditional knife-edge-Schliren in that the eye is much more sensitive to color changes 
than intensity variations, and there are now techniques to quantify the refractive index 
field. Rainbow Schliren permits a quantitative determination of temperature gradients in 
the gas phase, temperature gradients in the test section will deflect the beam and cause it 
to pass through different colors on the filter which will be recorded in the image on the 
camera. Depending on the orientation of the color filter, either horizontal or vertical 
gradients are measured. 

The basic equation governing the deflection of a light ray passing through a 


temperature gradient is: 


where dy is the deflection at the Schliren filter plane, f is the focal length of the 
decollimating mirror, L is the path length through the medium, On/OT is the change in 
index of refraction with temperature, with n, the reference index of refraction, and 
OT/Oy is the temperature gradient. Temperature gradients will result, therefore, in beam 


deflections that will be recorded by the camera. 


c. Heat Transfer to the Walls 

Standard Schmidt-Boelter type thermopiles together with surface thermocouple 
measurements have been used in most of the studies that deal with heat transfer at the 
walls and that are reviewed above [29-58]. The thermopiles provide linear millivolt 
ouput directly proportional to the heat flux through the gage. These gages have relatively 
fast response times (order .3 seconds) and are small in size (3.13 mm in diameter). They 


can operate in relatively hot environments (calibrated to over 473 K) and can be subjected 


to heat fluxes as large as 0.3 MW/m2. The response time and robustness of the device 


seem to be appropriate for the present application. 
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6. Characteristics of the Poitiers’ Intermediate-Scale Fixture 


A study on cross-flow effects over a pool fire has been recently developed at 
Poitiers and the main results are presented in references [71-76]. They mainly concerns 
have been the effects of confinement and cross flow on an intermediate scale pool-like 
diffusion flame. Especially the case where the fire is confined in a tunnel of similar 
characteristic length scale as that of the flame. The geometrical characteristics of the 
flame have been determined for different cross-flow velocities, fuel injection surface and 


velocity. 


(0.25 x 0.4 m’) 


Figure 2 — Experimental apparatus at the University of Poitiers 


A schematic of the diagram of the wind-tunnel is presented in figure 2. The 
experimental apparatus consist of a 5.25 m long horizontal channel through which air 
motion was induced by means of a centrifuge fan. The center part of the wind tunnel is 
the test section, 1.75 m long, 0.7 m height and 0.4 m wide (2 on figure 2). The walls, 
floor and ceiling are each made of seven removable actively cooled panels (0.7 m x 0.4m 
for the walls and 0.4 mx 0.25 m for floor and ceiling). The wall panels can be substituted 
by windows, for observation, and the floor panels by burners. The burners covers the 


entire cross section of the tunnel. Each burner consists of a rectangular sintered bronze 
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slab (20 mm in thickness) with imbedded copper cooling tube (8 mm in diameter). The 
thickness of the slab and the diameter are chosen to guarantee homogeneous flow through 
the entire burner surface. The burner width can be changed from as small as desired up to 
cover the entire channel width (0.4 m). The burner length can be varied in a discrete way, 
0.25 m, 0.50m, and 0.75 m, being the length the stream wise direction. Water from a 
reservoir at constant temperature of 65°C was pumped through the cooper tube to keep 
the burner surface at a given temperature and to avoid condensation. Similar burners have 
been successfully used in the past to simulate the pool or wall fire combustion of a liquid 
or of a solid fuels [71-73, 77, 78]. 

The characteristic velocity of the main flow can be obtained by means of a screw 
gauge placed in the symmetry plane of the test section or by Laser Doppler Velocimetry 
measurements. The wind tunnel can operate on both, blowing and aspiration modes. The 
rotating frequency of the fan is regulated to keep the flow velocity constant through out 
an experiment. In the aspiration configuration the air is conducted into the test section 
through a convergent (1 on figure 2), 1.8 m long, with a restriction ratio of 3.57. To 
homogenize the flow, a honeycomb, 50 mm in thickness, is fixed at the entry of the 
convergent. Combustion products are driven away through a collector (3 on figure 2), 1.7 
m long. To minimize environmental perturbations present in the laboratory, the floor of 
the channel is placed 0.9 m from the room floor. A similar arrangement is used when in 
the blowing mode. 

An initial aerodynamic study of the wind tunnel pointed out that the best method 
to induce a fully developed velocity profile upstream of the burner is by means of 
aspiration. Aspiration and blowing experiments were also conducted with flames and 
showed no significant differences. The geometry was chosen to guarantee, keeping 
minimum dimensions, that the flow will reproduce a turbulence intensity (u'/U) typical 
of fire scenario (between 5% and 15%) [72]. 

Temperature measurements can be performed using a combination of 50 um type 
K thermocouples and video recordings. A CCD camera records the visual characteristics 


of the flame on a U-matic video recorder. The image processing used to extract the 
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geometric parameters of the flame (length, height, angle of inclinaition, etc...) is a 
refinement of the technique extensively described in references [71-73]. 

Three parameters can easily be varied to performed a parametric study: the fuel, 
the cross flow velocity (U,,), the fuel injection velocity (V,), the the burner lenth (X,) and 
width (Y,). The common range for the air flow velocity (U..) is between 0.5 and 2.5 m/s. 
The fuel injection velocity is restricted by the source Froude number which, for fire 
flames, should be kept in the range of 10*-10°. The values used for the fuel injection are 
then between 1.9 and 5 mm/s. If the fuel is propane characteristic injection velocities will 
provide theoretical heat release rate is ranging between 18 kW and 135 kW. 

An example of the use of this experimental facility is presented in figure 3. The 
figure presents as presence probability distribution obtained from digitalization of a top 
view of the flame (figure 3(a)). From this image a mean flame is defined as where the 
presence probability is 50% (figure 3(b)). The zone where the intermittence is 50% is 
finite, therefore to obtain the main flame border the arithmetic average of the internal and 
external boundaries of the 50% intermittence zone is taken and define as the main flame 


border. 


(a) (b) 


Figure 3 (a) Presence probability of the flame in the transversal plane. (b) Mean flame contour. 
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It is not the intention of the present report to describe the experimental results 
obtained with this facility and that are described in detail in references [71-76]. The 
objective is to present the facility and its potential use for controlled intermediate scale 
validation of the small scale results derived from the bench scale screening test. 

Placing a simulated rectangular pool burner inside the horizontal channel of the 
wind tunnel describe above, it has been possible to efficiently manage air entrainment 
into the flame and to reproduce different levels of confinement. It has been underlined 
that the total mass flow generated by the flame is independent of the level of confinement 
and the flame height is significantly affected by a preferential orientation of the air 
entrained. Simple modifications to this facility will provide a “worst case” intermediate 
fire that follows well the philosophy driving the design of the bench scale screening test 


for SPGG’s, thus, provide a validation mechanism at a more realistic length scale. 


7. Numerical Methodology 


As mentioned before, the methodology proposed by McGrattan and al. [9] was 
used in the present program. The main characteristics of their approach to field modeling 
fire phenomena which emphasizes high spatial resolution and efficient flow solving 
techniques are presented in the above mentioned reference. This approach has been 
fostered by the ever-increasing power of computers and the development of faster 
numerical algorithms. To make the most of the current generation workstations, they 
have focused their efforts on developing relatively simple numerical algorithms to 
address the transport of combustion products in relatively simple enclosures. This 
approach will enable the use of a PC-Pentium computer for the computations. 

The mixing and transport of combustion products is calculated directly from an 
approximate form of the Navier-Stokes equations. This approximation involves the 
filtering out of acoustic waves while allowing for large variations in temperature and 
density [79]. This gives the equations a highly elliptic eheracter consistent with low 
speed, thermal convective process. In fact, a Poisson equation is solved at each time step 


in the calculation with an FFT-based direct solver [80]. 
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The fire itself is prescribe in a manner consistent with a mixture fraction based 
approach to combustion, but the combustion themselves are not, up to now, simulated 
directly and sub-grid scale motions are treated following the analysis of Smagorinsky 
[81]. The extent to which the fire is resolved depends on the objective of the calculation. 
As shown in the literature review, in a recirculating zone, such as the one of interest in 
this program,the reacting zone can be modeled simply as a source of heat and mass. The 
heat release and species production rates would be determined from experiment. 

Present capabilities permit fully three dimensional simulations at about 2 mm 
resolution for the proposed experimental facility. The products are simulated by tracking 
a large number of Lagrangian elements that originate reaction zone. These same 
elements carry the heat released by the fire, providing a self consistent description of the 
product transport at all resolvable length and time scales. Large temperature and pressure 
variations are permitted, subject to the limitation that the Mach number is much less than 
one. The methodology has already been successfully used to study a number of fire- 
related phenomena as shown in references [82], [83], and [84]. 

The LES code has already been used to perform successful preliminary isothermal 
flow calculations in two and three dimensional configurations similar to the proposed 
geometry (see Figure 1). I needs to be noted that the treatment given to the viscosity by 
the Smagorinsky [81] model results in an erroneous velocity determinations for a 
potential flow. In contrast, for a fully developed entry profile were viscous stresses are 
present the results are satisfactory. A simple correction for the zero gradient limit should 


be able to solve this problem. 
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8. Conclusions 


The literature review indicates clearly that is possible to well describe the 


structure of a recirculation zone located behind a backward facing step and consequently 


to design a fixture with which will allow to establish a diffusion flame stable over a large 


range of experimental conditions. Moreover, one can obtain the heat transfer profile at 


the burning surface and then, localize the hot spots were re-ignition would be favored. 


The following issues need to be taken into consideration when proceeding to the 


design of a test fixture around the recirculation zone induced by a backward facing step, 


baffle or trench: 


Mixing of fuel and oxidizer is guaranteed since splitting of the shear layer at 
reattachment results in part of the flow being deflected upstream into the re- 
circulating flow region to supply the oxidizer entrainment into the fuel zone. 

A significant increase in mixing occurs inside the recirculation zone since the flow on 
the entrainment side is highly turbulent, as opposed to the low turbulence-level 
downstream of the reattachment point. Care has to be taken not to increase the strain 
rate inside the recirculation zone to levels were extinction is induced by the main flow 
field. Systematic study of the operation conditions needs to be conducted as a 
function of the fuel to be used. 

The maximum measured back-flow velocity is usually over 20% of the free stream 
velocity. This value sets a criteria for fuel injection to the recirculation zone to 
guarantee optimal burning conditions. 

The flow in this region is very unsteady with very large turbulent structures (at least 
as large as the step height) passing through the reattachment region. This observation 
is of great importance in the present application since it ensures the entrainment of air 
to the re-circulation zone, adequate mixing and the evacuation of combustion 
products. The transient evolution of the recirculation | zone needs to be further 
evaluated for the conditions typical of the present application. The choice of a 


specific geometry such as backward facing step a baffle or a trench will have a 
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significant effect on the stability of the reaction zone. Stability of the recirculation 
zone results in low entrainment and low product evacuation, thus, not necessarily in a 
more stable reaction zone. 

e The aspect ratio of the test section has a negligible effect on the re-attachment length 
for aspect ratios greater than 10. For values smaller than 10 the reattachment length 
increases for a laminar boundary-layer and decreases for a turbulent flow. 

e Upstream perturbations can play an important role on the characteristics and stability 
of the eddy behind the step or obstacle. 

e Blowing or suction tunnels behave in a very similar way. 

e At least for high free stream velocities, injection has a de-stabilizing effect. No 
experimental data was found for velocities in the range of interest. 

e The length of the recirculation increases with Reynolds number, reaching a maximum 
at transition and then decay to a shorter length at the turbulent range. For Re>1500 
the recirculation zone has an almost constant length. For Re>6000 other recirculation 
patterns appear at the top and bottom of the channel increasing the complexity of the 


flow. 


Numerical calculations incorporating the k-e and algebraic-stress model over predict 
experimental observations and especially the turbulence dissipation rate. Preliminary 
computations conducted using the methodology described by McGrattan et al, on a one 
aspect ratio section, show qualitative agreement with the experimental data, further 
exploration of the use of this model needs to be done. 

For the low Reynolds number regime (Re,<2,000) flow visualization by means of 
illumination of particles introduced in the stream has proven to provide best results. 
Partivle Image Velocimetry is recommended as an appropriate technique to resolve the 
flow field characteristics. For punctual measurements of the turbulence intensity Laser 
Doppler Velocimetry will be necessary. Recommendations of seeding particle size and 
type are provided in the text of the report. Classic temperature and heat flux 


measurements seem adequate for the present objectives. 
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The test facility in Poitiers (France) seems well suited for future validation, at a 
more realistic length scale, of the screening procedure for SPGG’s. The information 
provided on the capabilities of this test facility can be used as background information on 
the design process. It also shows the potential of flow visualization techniques, velocity, 
turbulence and temperature measurements as tools to characterize air entrainment and 


flame structure. 
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